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The effect on flavour oscillations of simple expanding background space-times, motivated by some 
D-particle foam models, is calculated for a toy-model of bosons with flavour degrees of freedom. 
The presence of D-particle defects in the space-time, which can interact non trivially (via "particle 
capture" ) with flavoured particles in a flavour non-preserving way, generates mixing in the effective 
field theory of low-energy string excitations. Moreover, the recoil of the D-particle defect during 
the capture/scattering process implies Lorentz violation, which however may be averaged to zero 
in isotropic D-particle populations, but implies non-trivial effects in correlators. Both features 
imply that the flavoured mixed state sees a non-trivial flavour (Fock-space) vacuum of a type 
introduced earlier by Blasone and Vitiello in a generic context of theories with mixing. We discuss 
the orthogonality of the flavour vacua to the usual Fock vacua and the effect on flavour oscillations 
, in these backgrounds. Furthermore we analyse the equation of state of the Flavour vacuum, and flnd 

that, for slow expansion rates induced by D particle recoil, it is equivalent to that of a cosmological 
constant. Some estimates of these novel non-perturbative contributions to the vacuum energy are 
made. The contribution vanishes if the mass difference and the mixing angle of the flavoured states 
vanish. 
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I. INTRODUCTION AND MOTIVATION 



It has been suggested by Blasone and Vitiello and collaborators [1] that in quantum field theoretic systems with 
mixing, such as neutrinos with flavour mixing, for instance, the "flavour" eigenstates can be expressed as a condensate 
^ ^ in the Fock space of the mass quanta. Furthermore, it has been shown recently that there is a "vacuum energy" 
^ associated with this condensate which would be a contribution to the dark energy. In contrast to the standard mass- 
, eigenstate formalism, where the vacuum expectation value (vev) of the stress tensor vanishes, when one calculates the 
vev of the stress tensor with respect the flavour vacuum, the result is non-trivial, indicating spontaneous breaking of 
I , Lorentz-symmetry. In fact, computing such vevs for the Tqq and Ta components of the stress tensor in Robertson- 
^ ■ Walker space times, the authors of [2] have conjectured that there is a non-trivial, and non-perturbative, dark energy 
T-H ' contribution to the vacuum, and an equation of state p = wp with w < close to —1. These are features that are 
: compatible with recent observations. 

' In view of the breaking of Lorentz-symmetry by the flavour vacuum, characterized by a non-zero vev of the stress 
tensor, we will examine in this context a microscopic model discussed by Mavromatos and Sarkar [3] where such a 
feature is accommodated within a specific space-time foam model in the framework of modern brane/string theory. The 
foam model involves brane worlds which are embedded in higher-dimensional bulk space-times, and are punctured by 
DO-brane (point-like defects) [4] . The interaction of generic stringy matter with such defects involves the discontinuous 
process of splitting of the matter string by the defect, and subsequent joining, ignoring for the moment constraints 
from conservation of specific quantum numbers . This yields Lorentz-violating metric distortions [5, 6], g^^ oc u^u^, 
with the recoil velocity of the localised D-particle (DO) defect, due to the "topologically-non-trivial" interaction 
with the stringy matter as well as an induced local expansion. In a situation where one encounters a statistical 
population of DO-defects, one may have an average {{ui)) — (i.e. isotropic models for D-particle foam). Hence, since 
w^w^ = 1 for a velocity four-vector, one has a preferred time-like direction in space time in this case, leading to induced 
Lorentz-symmetry breaking from non-trivial higher order correlators of u^. Such a mechanism for Lorentz-symmetry 
violation is testable since it leads to a quite distinct signature for CPT violation [7] in meson factories [8], associated 
with a modification of the Einstein-Podolsky-Rosen(EPR)-type particle correlations. Electric charge conservation 
implies that charged states cannot undergo such interactions. Moreover "flavour" need not be conserved in the string 
splitting process. If D-particle foam contributes to the mixing of flavours then it only operates on the electrically- 
neutral neutrinos (or mesons). The above-mentioned Lorentz-violating "flavour- vacua" appear to be the physically 
correct ones to characterize the interaction of electrically neutral flavour matter (which for brevity we will refer to as 
neutrinos) with DO defects. 

It should be pointed out that in the work of [1, 2] the flavour vacuum formalism applies in general to all theories 
with mixing. In the above specific framework of stringy models of D-particle foam, only the neutral fiavour particles 
are subject to mixing and contribute to the dark energy. Dark energy, when interpreted as a cosmological constant, 
gives a de Sitter space-time with an associated horizon. Hence, it is interesting for self-consistency to re-examine 
and reinterpret the vacuum oscillation phenomena in simple space-times with horizons. In our case the latter are 
interpreted as the metric distortions due to D-particle recoil and will be shown to be most important for low energy 
neutrinos. This is to be contrasted with the approach of [2] where an a priori flat space-time formalism was used to 
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discuss the derivation of the dark energy contribution due to mixing. For simphcity we consider a two dimensional 
space-time metric as a physically-relevant toy model to illustrate the basic features of the approach. The restriction 
in space-time dimensionality is plausible if the interaction and recoil of neutrino matter with the DO brane defects 
(D-particles), is mainly along the (spatial) direction of the incident neutrino beam. 

The paper is organised as follows. We consider the recoil of D particles after the capture and subsequent emission 
of stringy matter in terms of suitable vertex operators [5, 6] within a logarithmic field theory [9]. By considering 
an ensemble of defects we model the metric induced by the recoil as a simple two-dimensional space-time with a 
horizon [10]. Finite systems on a lattice with periodic boundary conditions will be considered and the thermody- 
namic limit will be taken subsequently. This allows a rigorous treatment of the orthogonality ( and also of unitary 
inequivalence [1, 11]) of the flavour and mass vacua. The effect of the space-time expansion on neutrino oscillations is 
evaluated. In order to check if a dark energy interpretation of the equation of state is allowed, the equation of state 
of the flavour vacuum (after expansion) is calculated. Given the brane/string motivation for the model, a novel form 
of normal-ordering is introduced which demonstrates a contribution to dark energy. It is important to notice that the 
type of equation of state obtained by the above normal ordering procedure turns out to be formally independent of 
the initially assumed form of the scale factor of the Robertson- Walker background space time, thereby allowing for 
a self-consistent treatment of the back reaction non-perturbative effects of the flavour vacuum onto the space time. 
Some technical details of our approach are outlined in an Appendix. 



II. STRING INSPIRED TOY MODEL 



The discovery of new solitonic structures in superstring theory has dramatically changed the understanding of target 
space structure. These new non-perturbative objects are known as D-brancs and their inclusion leads to a scattering 
picture of space-time fluctuations. The study of D-brane dynamics has been made possible by Polchinski's realisation 
that such solitonic string backgrounds can be described in a conformally invariant way in terms of world sheets with 
boundaries [12]. On these boundaries Dirichlet boundary conditions for the collective target-space coordinates of the 
soliton are imposed. The properties of D-branes in curved backgrounds is not fully developed and so it is necessary 
to be guided by approximate calculations. In the interaction of low energy matter (e.g. a closed string) with a heavy 
D particle (i.e. embedded in a (c? -I- l)-dimensional space-time the D-particle recoils. Such recoil fluctuations of DO 
brancs can be described by an effective stochastic space-time. From our point of view there are two parts to this 
interaction process: precursor and aftermath of the capture of a string by a D-particle. Firstly the stringy matter 
has to approach the D-particle. In a superstring approach by considering the annulus amplitude for an open string 
fluctuation between a Dp and Dp' brane which have a non-zero relative velocity it can be shown that there is an 
attractive potential between them [4] . By specialising to p = 1 and p' = wc can give arguments as to why the low 
momenta strings are more likely to be captured by the D-particle. Secondly we have the capture and the re-emission 
of a string by the D-particle. This involves a change in the background of the string and as such breaks conformal 
invariance. It is not known how to describe the process of capture for short times. It is the capture and re-emission 
process which can be responsible for a change in the flavour of the matter string. For long times after the capture 
and emission process a semi-classical approach is expected to be valid [13]. In the wake of the capture process the 
vertex operator of the D-particle is described by an impulse approximation. The breaking of conformal invariance for 
large times can still be dealt with in terms of conformal data but with a set of recoil vertex operators which satisfy 
a logarithmic conformal algebra. Such an impulse/recoil event results in general in a time-dependent induced metric 
distortion and vacuum energy. It is such a metric distortion that is the motivation for considering the simple model 
that we will consider. 

The time scale of the D-particle matter interaction can be estimated from perturbative string theory. In the 
adiabatic approximation for the relative speed u between a p-brane and a D-particle separated by a distance r, it 
is possible to estimate the potential energy V{r,u) using perturbative string theory [4]. Typically the short range 
attractive part of the potential V is given by 

V(r,«) (2.1) 

with is being the string scale and is the only relevant piece in a superstring theory. For the case of flavour oscillations 
with sharp momentum k (whose associated mass eigcnstates have masses m^^^ and m^^^ with m^^^ > m^^^ in our 
convention ) the time scale r for the flavour changing capture and subsequent release of a matter particle by a heavy 
D-particle (when fi = c=l) can consequently be estimated by 
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where AV = 



_ 27r 



7^ 



7^ 



_ 2vk Am 



(1)2^(2)2 and Am^ = m^^^^ — m^^^^. Hence we deduce that quaUtatively 



4 m(i)^m(^)2 1 
27r Am2 A;^ ' 



(2.2) 



We shall then assume that the capture process is enhanced for low k. We can give an estimate for this time scale on 



considering a typical situation k 
is effectively instantaneous. 



IGeV, m(i) 



lO-^eV and Am^ ~ IQ-^eV^; we find t ~ IQ-'^Osec which 




FIG. 1: Schematic representation of the capture/recoil process of a string state by a D-particIe defect for closed (left) and open 
(right) string states, in the presence of D-brane world. The presence of a D-brano is essential due to gauge flux conservation, 
since an isolated D-particle cannot exist. The intermediate composite state at t = 0, which has a life time within the stringy 
uncertainty time interval St, of the order of the string length, and is described by world-sheet logarithmic conformal field theory, 
is responsible for the distortion of the surrounding space time during the scattering, and subsequently leads to Finsler-type [14] 
induced metrics (depending on both coordinates and momenta of the string state) and modified dispersion relations for the 
string propagation. 



We, for definiteness, consider a D3 brane world with D-particles in the bulk scattering off the brane. The recoil of 
a D-particle on scattering with a closed (or open) string (see fig. 1) leads to a world sheet deformation given by the 
semi-classical vertex operator V 

V = [ Gijy' {X°) dnX^ (2.3) 
■/as 

where Gij is given by the spatial part of the space-time metric, the world sheet boundary, 9„ is a derivative on 
the worldsheet normal to the boundary, is a spatial target space field satisfying Dirichlct boundary conditions 
and X° is the timelike target space field satisfying Neumann conditions (in the standard string a model) and y^{X^) 
is the classical D-particle trajectory which, with a suitable choice of co-ordinates, vanishes before the time to of the 
impulse. Let fci (^2) be the momentum of the propagating closed-string state before (after) the recoil, yi be the spatial 
collective coordinates of the D particle, gs (< 1) is the string coupling and Mg is the string scale and V^e~^{a' being 
the Regge slope) be identified with the target Minkowski time X^ for X^ after the collision For an initially 
Minkowski space-time background, we have [5, 6]: 

ViiX^z)) = Sijy\X\z)) = d {z) + Vi (z) (2.4) 

where z = e~*('^+*'^), (c, t) being the standard world sheet co-ordinates, and the operators Cj (z) and Di (z) are given 
(in terms of the regularised Heavyside function 6^ and Ui = jf- (fci — fc2)j ) by 

Ci{z) = syiee{x°{z)/V^) (2.5) 

and 

Vi {z) = UiX° {z) Qe {z) /V^) . (2.6) 

We define 0^ through 

( {z) /y/^) r -^e'o^°^^y^, e^0 + . (2.7) 

^ ^ 27^^ J_^ q-ie 
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We shall from now on choose units with a' = 1. These vertex operators have been calculated for non-relativistic 
branes where Ui is small. Consider now the operators C = CidnX'^, and V = VidnX^. By calculating the correlators 

{C{z)Ciw))^0{e'), (2.8) 
(C(z)PH>~^-^, (2.9) 



\z — w\ 



and 



L 

\z-w\' 



{V{z)V{w))^^-^^-^ (2.10) 



with respect to the measure of the free bosonic string action we see that we have a logarithmic conformal field theory 
[9] , provided that [5] ^ (log ^) , where L is the diameter of the world sheet disc and a is an ultraviolet world-sheet 
cut-off. Hence e is non-zero for finite L. 

By considering the operator product of these operators with the stress energy tensor we find an anomalous dimension 
A = ^ . In order to restore conformal invariance we need Liouville dressing with a dilaton field. For simplicity if we 
restrict the slow recoil to the x direction (on the brane) we obtain the vertex operator V^'' 

V^^) = / d^z e"°«^a„ {u^X^ {z) Oe (X° {z)) {z)) (2.11) 

is 

where (j) is the dilaton field. From considerations associated with conformal invariance of the bulk world sheet of the 
dressed theory [15] 



ao. = -| + ^^ + y (2.12) 

where Q is the central charge deficit of the cr-model describing the stringy excitations of the recoiling D-particle. We 
remind the reader that the restoration of the conformal invariance by the Liouville (f) system implies 

c^ = 1 + 3Q2. (2.13) 

where is the central charge of the Liouville part of the a- model. As discussed in detail in [5, 6], one may argue 
that, in our case, is O (e'*), assuming a critical bulk space-time dimension (c.f. fig. 1); hence, to leading order in 
e 0+, its contributions can be ignored in comparison to the anomalous dimension term in (2.12). Thus, we 
have aox ~ (for u<^l) and from (2.11) we obtain: 

V^"' « / eTf (z) {z) u^Oe {x" {z)) + ... (2.14) 

is 

where the . . . denote terms that do not contribute to physical correlators. Such terms are proportional to either 
X^S{X^) (obtained by the action of da on the Heavyside operator Qe{X^)), or dad°'X^ (which vanish in the conformal 

gauge on the world sheet). 

In the recoil case the string is slightly supercritical, i.e. > 0. In such a case, the zero mode 0(o) of <p can be 

interpreted as a second time in addition to X° [16]. From (2.14) we then observe that the effect of V^^^ in the action 
can be reinterpreted as an induced target space metric in the extended D + 1-dimensional space-time (^(0),-^°,-'^') 

€4,4, = 1, C = -1, Goo = -1, Go. = eV5 (X") (2.15) 

where, following [13], we have assumed euclideanised as required for convergence of the world-sheet path integral. 

At first sight, it appears that this metric is characterised by two times, ^(o) and X^. However, it has been argued [13] 
that the open string excitations on the brane world are such that they live on a D-dimensional hyper-surface of the 
extended D + 1-dimensional two-time space-time; the Liouville mode is proportional to y^^,-o M^^^,with the 

four- velocity on the D3 brane for the recoiling D-particle and can be written as 7(1, w,), where 7 (= l/^/l — 
is the Lorentz factor . 
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Wc shall demonstrate this structure in an explicit example in which the threc-brane world is obtained by compact- 
ification on magnetized tori, characterised by a magnetic field with intensity H. This magnetic field is not a real one, 
but a "statistical" field characterising the compactification. An open string excitation coupled to such a magnetic 
field, has its ends attached to the compactified three brane world. An encounter of the string with a D-particle is 
described by the addition of the following deformation to a free open-string cr-model action: 



HQ^{u^X'')X''d^X'^ (2.16) 

as 

where denotes tangential derivative on the world-sheet boundary. This corresponds to the vertex operator for 
a gauge potential corresponding to a uniform magnetic field in the compactified toroidal dimensions {X'^,X^), that 
appears suddenly at a time X^^^g '^ij-^^ = 0) co-moving frame of the recoiling D-particle, which has captured 

the open string excitation. The reader should notice that the specific calculation takes place in this co-moving frame; 
the argument of the Heavyside operator is the rest- frame time of the recoiling massive D-particle, which in terms of 
the original coordinates is expressed as u^X^^ = j{X'^ + WjX*). 

It can be easily seen [5] that the operator (2.16) has anomalous dimension — ^u^u'^ ~ ~T' since the four velocity 
satisfies by definition u^u'' = 1, with respect to the background Minkowski metric. This implies a Liouville dressing 
similar to (2.11), with the corresponding Liouville anomalous dimension a 

[ e^'^He,{u^X^')X'^^''X\ (2.17) 
Jdi: 

As shown in [13] the open superstring excitations that couple to such a magnetic field will exhibit Zeeman-like 
supersymmetric mass splittings between fermions and bosons of the form: 

dm^ ~ i/eTf '^iJee(w^X^)E45 (2.18) 

with E45 the spin operator given by j [74, 75] where 7^ (/z = 0, . . . , 5) are the Dirac matrices relevant for 6 dimensional 
space-time. However masses need to be time (or equivalently e) independent for stability. Representing the Heavyside 
operator for positive arguments as Q£{u^X^) e"^""'^ , Sm^ is independent of the value of s if [13] 



7f " 



Uf,X^ = constant (2.19) 



where the constant can be taken to zero without loss of generality. Since the above-mentioned mass splittings can 
be tuned by adjusting appropriately the intensity of the "magnetic" field H, in a way independent of the size of 
the four-dimensional vacuum energy contribution [13], the above example is sufficiently generic and so of interest for 
phenomenological applications of superstring theories. 

The constraint (2.19), which is consistent with general coordinate invariance on the brane world, implies that the 
dynamics of the open string excitations can effectively be considered on a D-dimensional space-time (hypcrsurface), 
with a single time. Although above we have argued in favour of (2.19) in the context of a specific superstring 
situation, nevertheless we may assume its validity in more general frameworks of supercritical strings [17]. 

Prom (2.19) and (2.15) we then obtain for the induced Z)-dimensional target space-time metric (with X° = t): 

ds^ = #2 _ dt^ + u^dXdt - {dXf = 27^(^4 -I- u^dXf - dt^ + u^dXdt - {dXf (2.20) 

We next consider the case of a statistically significant population of D-particles on the brane world (c.f. fig. 1) 
("D-particle foam"). We may assume a random (Gaussian) D-particle foam, such that [3]: 

((«,)) = , {{u^u^)) = a^, (2.21) 

and we assume slowly moving D-particles, so that cubic or higher powers of the recoil velocities are ignored. The 
average of the metric (2.20), then, yields (upon expanding appropriately the Lorentz factors in powers of u^): 

((ds^)) = (1 + 2a^) dt^ - (1 - 2a'^)dX^ + dY^ + dZ'^ (2.22) 

Notice the existence of horizons in this metric, the "xx" metric component vanishes for = 1/2 (which is consistent 
with the fact that the velocities are subluminal \u^\ < 1 ). The metric (2.22) has a global character, as a result of 
averaging over populations of D-particles that cross the entire D-brane world. 
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We next notice that the quantities are in principle time dependent, denoted by cr^(t), since the density of 
the D-particle defects may depend on time. This can happen, for instance, if the bulk density of the D-particles 
is not uniform. The nature of the time dependence of (J^(t) is therefore dependent on a microscopic model. If the 
density of the D-particles that cross our brane world (c.f. fig. 1) decreases with time, we then have an expanding 
Robertson- Walker type Universe in the directions, and a static universe in the transverse directions ^.Z): 

{{ds^)) = _ a^i^^)(^dX)^ + ... (2.23) 

where the "cosmic time" <r is defined as c?<^^ = {l + 2a'^{t))dt'^ and the . . . denote the static directions which will not be 
of relevance to us here. It is a priori natural to consider an isotropic D-particle foam, in which the D-particles recoil in 
random directions so that {{ui)) = 0, {{uiUj)) — cr^Sij . Asymptotically a(?) monotonically tends to a constant with <;. 
In such a case one obtains an isotropic Robertson- Walker like Universe. For simplicity in this work we concentrate on 
effective two-dimensional universes, since such an assumption will not affect the qualitative features of our approach. 
In this sense, the interaction of D-particle foam with stringy matter on brane worlds, leads to an effective space-time 
metric, "felt" by the matter, which, depending on properties of the foam, can have the form of homogeneous and 
isotropic Robertson Walker cosmology. 

III. EFFECTIVE TOY MODEL FIELD THEORY OF BOSONS WITH MIXING 

We will use the string theory considerations, leading to (2.23), as a motivation for studying the possible consequences 
of our D-particle/string hypothesis for flavour mixing within the context of a simple field theoretic 1-1-1 dimensional 
model with a horizon. Before starting, we mention that flavour need not be conserved during the capture process of 
fig. 1, in the sense that the re-emitted string, after the capture, may have a different flavour (and mass) from the 
incident one. However, the spatial momentum is conserved in the process. In this sense, our D-particle foam may be 
considered as a medium of inducing flavour mixing [3, 18], which might be a feature of quantum gravity foam^ in 
general. 

As we shall discuss below, the effects of flavour changing on the space time arc drastic, and resiilt in a "new Fock- 
space vacuum" state of the form suggested in [1] for generic theories with mixing. As we shall discuss in section V, 
this results in non-trivial contributions to the vacuum energy, of cosmological constant type, which come mostly from 
the infrared modes of the flavoured particles, of momenta less than the a representative scale of order of the sum 
of the masses. These constant (in time) contributions imply that, eventually, the flavour-changing process dominate 
the simple Robertson- Walker metric structure (2.22,2.23) arising from the D-particle recoil distortions, results in de 
Sitter type Universes. 

A. Formalism in an Expanding Universe 



Let us consider the metric in 1 -|- 1 dimensions [10] 

ds^ = dt'^ - c? {t) dx^ (3.1) 

which in conformal time co-ordinates has the form we have 

ds^ = C (ry) {drf ~ dx^) . (3.2) 

According to our discussion above, we first consider Robertson- Walker expanding Universes, interpolating smoothly 
between asymptotically flat space times. In our D-particle-foam picture, leading to (2.22,2.23) this can be achieved by 
considering models of foam in which the density of the D-particles in the bulk is not uniform but varies appropriately 
in such a way that, as the D3-brane world moves through the population of D-particles (c.f. fig. 1), the effective density 
of the D-particles crossing the D3 brane varies smoothly with time so as to lead to the above situation of smoothly 
expanding isotropic Universes. Moreover, we assume for simplicity the case in which the fluctuation parameter cr^ {t) 
never reaches the horizon scale 1/2, but is always kept sufiiciently small. This is only for simplicity, since the essential 
features of our expanding Universe, that would be relevant for our discussion in this work, namely particle production, 
can already be captured in such horizon-free situation. In any case, as we have already mentioned, and shall discus in 
section V, non-perturbative effects of our flavour vacuum will produce de Sitter space times with asymptotic future 
horizons. 
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A simple choice for C {rf), capturing the above features is 

C{ri)=A + B tanh [prf) (3.3) 

with ^ > S > and p > and has the feature of C (77) asymptotically tending to a constant as — > c». This space 
time is asymptotically accelerating, for sufficiently negative prj —00, and decelerating for positive prj ^ 00; hence 
there are no asymptotic future cosmic horizons. This is consistent with the absence of such horizons in perturbative 
string theory [19], employed in our world-sheet approach to D-particle foam, which is based on well-defined Scattering- 
matrix elements and asymptotic states. In the D-particle picture given above, if we assume that asymptotically in 
cosmic conformal time ry ^ 00, the recoil fluctuations tend to a constant value, ct^, we can identify A ~ 1 — 4a^ and 
B = 2a^. Of course, there is also a flat space limit when is B ^ 0, i.e. vanishing D-particle recoil. The space-time 
(3.3) is considered as a perf,urbatwe stringy background on which we shall discuss flavour changing process resulting 
from the stringy, topologically non-trivial processes of capture in our foam model, depicted in fig. 1. As mentioned 
above, siich fiavour changing process will lead to a new vacuum state, which in turn will imply contributions to 
the vacuum energy of the Universe, of de Sitter type, that will eventually dominate the space time (3.3), leading to 
eternal acceleration and asymptotic horizons. This should be considered as a non-perturbative contribution of the 
flavour-changing non-trivial vacuum structure of our model. 

On the other hand, if one would like to use an exactly soluble model, without future horizons, interpolating between 
a de Sitter space time at a certain era in the past and a flat Minkowski space time at the future asymptotic end of 
cosmic time, then (s)he could consider as an instructive example the following form of the scale factor [20]: 

C{r])=A + B^^= (3.4) 



with the model incorporating the de Sitter phase in the far past, pry ^ — 1 for A = J5 > 0, p > 0. Indeed, for 
sufficiently negative conformal times, such that (pj?)^ ^ 1, the model implies infiationary evolution, since in that 
range of parameters 

On the other hand, in the asymptotic future r]p ^ 1, the model asymptotes to flat Minkowski space time, with 
C(t]) 2A. 

In what follows we shall not discuss explicitly the space time (3.4) further, but concentrate instead on (3.3), since 
for our qualitative purposes in this section of demonstrating the eff'ect of expansion insofar as particle production is 
concerned, the two space times (3.4) and (3.3) are equivalent. However, we note that, on anticipating a small Hubble 
parameter of our D-foam-induced-de Sitter space time (c.f. (5.34)), we are essentially interested in the small p regime 
in which both space-times (3.3) and (3.4) undergo linear expansion with the (conformal) time. This limit of slow 
expansion, p — !■ will be understood in what follows. 

It goes without saying that the above issues are far from being considered as resolved, especially in a string 
theory context, where the issues of an inflationary phase, and in particular a smooth ("graceful") exit from it, are 
far from being understood analytically. In critical string theory, there are recent attempts towards a construction 
of exact solutions, incorporating inflation, by making use of the properties of the compact Calabi-Yau spaces in 
concrete, phenomenologically semi-realistic (as far as incorporation of the Standard model is concerned) brane world 
models [21]. However, the exit phase from the de Sitter era, together with attaining the very small value for the 
vacuum energy observed today, are still not understood, in our opinion. Some attempts towards these latter issues 
do exist, however, within the context of non-critical (Liouville) string framework, which the present D-particle foam 
model belongs to. For further details we refer the reader to the literature [22]. 

In the present article, we shall certainly not attempt to offer any quantitative resolution of these problems, given 
that our considerations below will be based on effective field theories. Nevertheless, we think that our findings are 
sufficiently interesting to probe further studies along this direction, which might yield useful contributions to the 
quest for a proper understanding of quantum gravity and its connection with other unresolved issues in fundamental 
physics, such as the fiavour mixing problem. 

After this necessary explanatory remarks, we now proceed with our analysis. In conformal co-ordinates our 
expanding-Universe background space time is written as 
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with |detg^i/| = g = C"^ (77). Wc will consider the dynamics of a real scalar free field ^ (appropriate for a neutral 
meson) of mass m in this background: 

C-^{v)d^{C{v)g''''d,<P) + m^cj> = (3.6) 

This simplifies to 

{d^^-dl)<j> + m^C{v)<i> = 0. (3.7) 
On writing [10] cp (jj, x) = U (??) V (x) we can find a complete basis set for expanding fields from 

^V + eV = (3.8) 

and 

-^U+{e+m^C{r]))U = 0. (3.9) 

We have V {x) = e*'^^ and U{r]) = ciUi{r)) +C2U2{r]) where Ui{r]) and U2{r]) are given in terms of the hypergeometric 
function F (a, b; c, z) (= Y^)r(b) ^n=o ^^"r(c+i^)"'""^ wf) • ^^^^ ^® helpful to introduce the notation 



CU^n,k = ^k^ + w? {A-B), (3.10) 
^out,k = ^/WT^^^?{A + B), (3.11) 



and 



'^±,fe = 2 {^out,k ± Wi„,fe) . (3-12) 
In terms of these variables it can be shown that [10] 

C/i (ry) (3.13) 



oc (1 +tanh(p?7))"'"' (1 - tanh (p??))''^" = F 1 + -a;+; 1 + i— , - (1 + tanh(p?7)) (3.14) 

VP P P 2 / 

where ai = ^ ^1 + ^^^^ ^ /J;^ = i + ^ . We can find another independent solution U2 {if) on making the change 
•q —rj and B —B (i.e. ujm < — > Wout) and so 

U2 (ti) (3.15) 

oc (1 - tanh(pry))"^-^ (1 + tanh(pr;))'^^-^ F ( 1 + 1 + i — , - (1 - tanh(pr?))^ . (3.16) 

\P P P 2 J 

We should note that complex conjugation of Ui (77) is also a solution of (3.9) and this is effected by taking the complex 
conjugations of the arguments of F and also the exponents in (3.13) and (3.15). In terms of these solutions the 
appropriate in and out solutions can be taken to be 

<(7?,x)=e*'=-C/r(r?) (3.17) 

and 

uT'{ri,x) = e'^-U2{ri) (3.18) 

where the constants of proportionality in (3.13) and (3.15 are chosen on the grounds of normalization. It is known 
that there are Bogolubov coefficients ak and (3k such that [10] 

4" (r?, x) = akur' (r?, x) + pkulT iv, x) (3.19) 
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and 



with 



and 



oik 



-r(i-^)r(- 



r 1 - 



^out,k 



.,k r(^)r(i 



(3.20) 



(3.21) 



(3.22) 



(N.B. the Umit of no expansion is given by B ^ with ak ^ I and — > as required.) This formaUsm requires a 
shght generahzation when there are 2 distinct mass eigenstates. We will need in the obvious way will, ^out^ with 
i = 1,2 with m replaced by m^*). In terms of creation {a^} and annihilation {a} operators for massive particles in 
the in and out Hilbert spaces we have the relations 



in(i) _ (i)* ouL{i) _ ouL 

out(i) (i) in (i) , ^(i)* m (i)t 

% =^k % +Pk «-fe • 



where for i = 1,2 



j^) r(i 

^out.k 



VS. r(-=4i)r(i-=£ 



(3.23) 
(3.24) 



(3.25) 



and 



r(i 

^out,k 



„(■) 



r 1 + 



The corresponding massive fields have the standard non-zero commutator 

(ri, x) , dt 4>^^^ iji, y) ^iS{x- y) Si. 



(3.26) 



(3.27) 



It will be useful when discussing issues of orthogonality for condensate states that we put our spatial system in a box 
of length L with periodic boundary conditions i.e. with k = ^jf-; in natural units 



</)W(r7,x)= ^ 



2Luj 



(3.28) 



with Lu'f^ 



(2|i) +m«'and so 

0W {rj, x) , dt<f>^'^ {rj, x')] = i le^^ = iA^ {x - x') . 



(3.29) 



We also introduce a short distance cut-off Sx (or lattice spacing) so that wc have a discrete lattice structure: NSx = -j 
where A'' is an integer. There would then be a maximum \j\ = N [— and so —N < j < N . Similarly with 
momentum —-^ and so — A'' < j < N. Space points become discrete: x nSx and x' n'6x so that — A'' <n,n'< N. 

We can then define A^^^^ [x — x') = Yl!j=-N e*^^"'^ )^ as a lattice analogue of {x — x'). Consequently on the 
lattice 



N 

{n6x) = (a; 



(3.30) 



with -N <n< N. 
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B. Flavour mixing and Vacuum 



Although it would be desirable to rigorously demonstrate string matter flavour changes due to D-particle capture 
this requires non-perturbative string field theory or M theory. Since such a theory is not available we will consider a 
field theoretic treatment where the mixing angle 9 is put in by hand. It is later when we consider the vacuum energy 
contribution due to D particles that we will require a correlation between 6 and D-particle recoil. Following [1] flavour 
mixings are generated by G (O) 



g [6) = exp {ies) 



where <S is given by 



3=0 



where 



(3.31) 



(3.32) 



(3.33) 



with 



7j± = 



,(1) 



± 



,(2) 



1) 



(3.34) 



and 



a'-'^e "^j *. There are of course additional superscripts associated with the in and 



(and so j-j± = "fj±, 
out Hilbert spaces. Clearly we have the standard commutation relations 



~t(i') 



~t(i) ~t(i') 



= . Hence 



[S,-,S,.]=0, 



6jj' dii' ,and 



(3.35) 



j=i V / 



and 



JV 



3 = 1 



It is then straightforward to check that 



g-^ (0) age (9) = cos^sg + ^ sine (• 



7i- '4^^ + lj+ 



and so 



5g) 



9-' (9) «gS (») = co.«a« + i=i„e (7,-a«f e*W"-^"3)* + 7,+«ge'W" -?')') . 
Similarly we have (with j > 0) : 

g-' (e)Sg^ W = cos^Sg + i sine (7,-4? - • 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



(3.40) 
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These relations for creation and annihilation operators are consistent with the following mixing relations for the 
corresponding massive fields (/>^'-*: 



g-^ (6) (r/, x) g (0) = cos e (j)^^^ {-q, x) + sin 6 (f>^^^ {r], x) (3.41) 
g-'^ (6*) {n, x) g {0) = cos ^^^^ (r?, x) - sm {rj, x) (3.42) 

the right hand sides of which can be identified with (rj, x) and $/3 (?7, x) fields with flavour quantum numbers a and 
/? [IJ.Thc annihilation operators „ and aj^p for particles with flavour a and f3 can bo identified with g-^ (9) a''pg {6) 
and g~^ {6) af^g {0) respectively. Since we are interested in the particular contributions of D-particles to mixing the 
above definition of flavour vacua |0)^ ^ will be restricted to k < k — ^'^^ 



a,- „ |0)„,^, = 0, (3.43) 
|O>„,0 = (3.44) 

on assuning that D-particle mechanism is responsible for particle mixing. For values of k higher than the cut-off |0)^ ^ 

will coincide with the massive vacuum i.e. a^!^- |0)^ ^ = for i = 1,2. This is an important difference between our 
model and that given in [1] . Hence we can deduce that 



g. (0) |0)„ . = |0), 2 (3.45) 



AT* 

where ^* {0) = JJ exp (— fSj) and so 



|0)a,/3=^*"'W|0)i,2- (3-46) 



For finite N* the flavour and massive particle vacua have a non-zero overlap. In the thermodynamic limit N* — > oo 
this ceases to be the case as will be discussed in the next section. 



C. Relations between flavour and mciss-eigenstate vacua 

We will be interested in the N dependence of the overlap / (0) between \0)i 2 ^^'^ l^)a,/3 

f{0)= i,2{^K'm\0)^,2 (3-47) 

which has a factorized structure since g~^ {0) = JJ^/^ where gj^ (0) = exp(|Sj). Hence it is sufiicient to 
consider 



fj{e) = ,,2{o\giHo) |o)i,2. (3.48) 



Now 



^/iW = ^ 1,2(0 |6,g7i(^) |0)i,2 



7i- 



(a«Sg+agsf)a7^(^)|0),,, (3.49) 



and 



{s^^^+^^)g-^i0) (3.50) 
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Also, since [&j,Qj'^ (0)] = 0, 



^fA&)-liAo\g-'{0)e,\o),^ 



i,.(0|e-^W(a«yf +sL^^f ^) 10),,,. (3.51) 



Prom 3.51 and 3.49 we deduce that 



1 + 2^ sin2 e] ^ fj {9)=-^ sin [26] f, (0) . (3.52) 



which has a solution 



„2 



/,((^)(x |^l + ^sin2(0)j . (3.53) 

Because (m^^^^ — m^^^^) is small, within the context of our model we can still choose a sufficiently large j (> A''' say) 
such that 

and so 

j=o j=JV'+i 1 + ^^'iyi ' sm^ j=o 1 + ^ sm 



where A^' ~ "^'^y^ 



AT' 

The factor — = Fjv satisfies the bound 



i^iv < n ^ TT (3-56) 



=il + a 



1 - ^ 



where a = j sin 9 - — — [sr^- The other factor in (3.52) can be written as 



n TTX 

where o = j-^ . Now we can define 



3 

and so 

JV* 



f)(^«) = n^^-7Y (3-57) 



l_ _ I) (5, A-*) 
J 

j=JV' + l " j 

As A'* — > oo also A'' — > oo and so we examine the asymptotic behaviour of f) (6, n). It can be shown that 



r2 (n + 1) r f 1 - r f 1 + i\fb) 

f) (6, n) = — . / V (3.59) 

r f 1 - iV6 + n j r n + iV6 + 



13 



Prom Stirling's formula as n — > oo 

r2 (n+l) 



(l-iVb- 



+ n]T 1+i 



i\fh + 



1 - 



1. 



JV* 



Thus JJ contribute to any orthogonality of the flavour and mass vacua. An upper bound estimate 

j=N' + l ^ 

^N' for Fn' given in (3.56), is 



5;v^ = (l + ar- 11—^ 

j=l l+a N' 



Now it can be shown that 

N' 

logH 



1 _ a 3 

j=l ^ l+a N' 2 



2 + log (1 + a) - 2 



a+l 



tanh ^ 



l + a 



and since a can reasonably be assumed to be very small and so 

JV' 



< N' 



1 7 

r-30 



a2 + o(ai) 



From (3.60) 



\^N' I < exp 



2a N'' 



(3.60) 



(3.61) 



(3.62) 



(3.63) 



and so ^Jat' ^ as TV' ^ oo. This demonstrates that the flavour and mass ground states are orthogonal separately in 
the asymptotic in and out Hilbert spaces provided a is non-zero ( i.e. for non-zero sin^ 6 and (m^^^ — m^^^ ) ). This 

o 2,2 

orthogonality is actually robust to modifications of the energy-momentum relation of the form = — ^ 7™ \n due 
to possible quantum gravity effects where n is a positive integer [23] . 



IV. VACUA AND THE EFFECT OF EXPANSION 



A. Orthogonality 

The effect of expansion is well known to result in particle production for co-moving observers. Consequently we 

shall consider its effect on the orthogonality of the flavour and massive vacua. The relation between in and out 
operators (c.f. (3.23) and (3.24)) connected by Bogoliubov coefficients can be effected in terms of an operator B. In 
terms of 



exp 



(i) (i)t (i)t (i)* (i) (i) 



and 



[f, J = exp y-i^) 'ay ^ «-V«-iJ j 



B i^^j' 1 ^^j*'^ can be written as S (^'-'j*'') P (''^■'*') ' shown 



that 



{%)in 



ii^out (i)out\ i%)out Tt ( (i)out (i)out\ 



{i)out 



■ sinh ( 



{i)out 



i(i)out , -, / (i)out \ 

a_j + cosh [ Vj j • 



3 a- 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



on identifying 



V. 



■{i)out 



e = a 



(i)* 



cosh (I ^^f""* I) 
where = v^^^"^*e^'^i . Analogously 

(i)oui 7->t / (i)in (i)in\ (i)in r-> ( (i)in (i) 

a)> = B^[ipy ,vy )ay B[^)> ,vy 



with 







■~^(i)out 









■ tanh ^ 



sinh 



sinh ( - 



.,<■) 



where we have used the natural generalisations of (3.10,3.11 and 3.12) defined by 



in,j 



47^2 j 2 



+ mW2 {A-B), 



and 



out J 



UJ. 



'±,3 



■mW2 {A + B), 



(i) 

out,j 



in,j J 



Also 



. (i)out 



out, J 



r 1 



J r ( ^ r 1 + ^ I cosh I 



/ ~{i)in \ 



Moreover v^^^ = — and ip^^" = — i^j*^""* follows straightforwardly. 



We have now the basic operators to calculate the overlap amplitude 
which is equivalent to 

fm(l,2;a,/3)= i,2,in(0|C'W|0U,i,2- 
The relation between in and out massive vacuum states (|0) -^^ ^ 10)^^^^ 2) can be deduced from 

since we require 
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Hence 



Similarly 



i=lj>0 



I out,l,2 ■ 



i=lj>0 



2 • 



(4.17) 



(4.18) 



i=lj>0 



where G^out (^) defined analogously as before (c.f.. (3.31)). 



l,2,out 



(01 n^( 



(i)out (i)out 



in terms of which 



fm(l,2;a,/3) = n^fcW- 



fe>0 



In the small approximation ( appropriate for the size of effects generated by quantum gravity) 



In terms of the operators 



V2 



V2 



(with j > 0) 



c/ (»)!,Wt2 0)*,W2\1 r 1/ (i),(i)t2 (i)*,(i)2\ 



(4.19) 
(4.20) 

(4.21) 

(4.22) 
(4.23) 

(4.24) 



which has the form of a product of standard squeezing creation operators S (^) = exp [| (^a^^ — ^*a^)] which generate 
a Fock space representation of squeezed states 



n=0 



where 



C=^tanh(|^|), 



On using this property, in terms of Cfe^^°"* = j tanh ^ ^ and C,j^ 

shown that 



(2)out 



(4.25) 
(4.26) 

■ tanh ^ y^^'>°^* ^ ^ it can be 



= \ 1 



j.(l)out 



J2)out 



(4.27) 
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' out ,j 



|0) 



out, 1,2 



0. 



Let us now compute the remaining terms in (4.21). We can show that 

i,2,o„t( 0| (t;f°''*)p(0f™*) 6, 

U=i J 

The second order term in 9 gives 

i,2,o„*( 0| P (^Wout) g2^^ .| |o^^^^^^^^ 



(4.28) 



-27 



1- c 



(l)out 



1 - 



(2) out 



and so 



•^i(^)^\/l- 



(l)otit 



1 - 



(2)out 



1 2 



1 _ _f>-2jr/o«t,3 
4 



^^(l)oMt* ^(2)out*y 



(4.29) 



(4.30) 



If the small 9 expansion had not been made the total expression would have been periodic in 9 and so, to be compatible 
with the earlier calculation ( in the absence of expansion), we can rewrite (4.30) as 



•^i (0) 



1 - 



(l)otit 



1 - 



(2)out 



1-1-^2 sin^e ^ _ 3^-2iT,out,3 ^^(l)o«**^(2)o«t' 



(4.31) 



which is again valid for small 9. Now the denominator can be bounded 



1+7, 



sm 



1 - Je— - (C 



r 



sm 



16 



4-3 



(l)ouf 



Hence 



\^A0)\'< 



1 - 



(l)out 



1 - 



(2)out 



(2)out 



1 + 7 



2 sin^ 9 
16 



4-3 



.(2)ou/, 



(4.32) 



For large j, 



^{i)out 



-e 



27rip 



which is a very sharp fall off with j. To leading order in B for small j and for small p ( i.e. slow expansion) 



exp 



2pAV2 



(4.33) 



(4.34) 



This behaviour for 



(i) 



does not change the previous arguments for orthogonality and so the flavour vacuum (after 

expansion) and the original massive vacuum in the thermodynamic limit are orthogonal. In [1] it was asserted that 

this orthogonality implies that the two vacua are unitarily inequivalent, i.e. there is no unitary transformation that 
connects them, and that the Fock-space states constructed from one vacuum are orthogonal to the Fock-space states 
constructed from the other. The lack of a unitary transformation connecting the two vacua can be explicitly seen [11] 
in our approach by observing that the flavour vacuum is not normalized to 1, as should be the case if this vacuum 
were unitarily equivalent to the standard mass eigenstate vacuum. 
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B. Oscillation Probability 



Particle production during the expansion will lead to a modified oscillation and we shall examine the signature of 

this modification within this model. In the presence of expansion the transition amplitude Oytoscio^jk'; (3,k) for the 
oscillation of fiavour a (momentum k' > 0) to /3 (momentum fc > 0) is 



Mosc [a, k'; (3, k, rf) = k'jni "I P iv)) 



out.k ' 



(4.35) 



where (by definition) \l (j?))o„j = af™* |0)otit a fi ^'^'^ i = a, (} and similar definitions apply for in states. (Previous 
formulae with discrete momentum labels j, when required, translate readily to the continuous k label.) Hence, in 
terms of the massive Fock space, 



Wosc {a, k'; 13, k, rj) = i,2,m( 0| jcos {9) a^.!'*" + (7ft' -a-fe'*" + 7fe%+ai^ 
^ / ^~(2)o«tt sin6» / 



-k - 7fe,+afe 



/out, 1,2 ■ 



It can be shown that 

|0) out, 1,2 
f 

= n 

fe>0 



1 - 



{l)in 



1 - 



»{2)in 



V 



E„+,,.„ Sr (2n+,i, 2n+,2, 2n_,i, 2n_,2) 
|2n+,i, 2n+,2, 2n__i, 2n_,2),„ 



where 



0*fc"(2n+,i,2n+,2,2n_,i,2n^,2) 



exp (21 liy^fe-*'" [n+, 



(2) 



}) 



^ ^(l)m [n+,i+r»_,i]^(2)m [n+,2+n-,2] V (2»+,l)! (2rt+,2)! (2n_,i)! (2?^^^ 



2n+,i+n+,2+n-,i+n_,2 (n+_i! n+,2! n_,i! n_,2!) " 



Hence 



fe">0 



.(1) 

Sfe" 



1 - 



a(2) 
Sfe" 



X <5fe'fc {x:i,kVi,k iv) + >i2,kV2,k (7?) + ^3,kV3,k iv) + >t:A,kVi,k {n)} 
with xi,fe = ^ cos^ = -X2,fe, X3,fe = 2^ sin^ 7fe,+ = -X4,fe and 

= ^ sine (7fe,-/?i'^*e-'4^'°"*'/ - 7fe,+a«*e'4^""'*'/) = (,?) 

and 

-U4,fc (77) . 



f 3,fe (?y) = ^ cos 6* ^ e'^fc 
v2 



(4.40) simplifies to 



Ttosc {a,k';0,k,ri) 



= ( n Jfi-M" 



a(2) 
''fc" 



X sin 20 



(4.36) 
(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 
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and can be readily interpreted as a modification of the usual formula for oscillations with the extra momentum 
dependence (due to the field theoretic treatment); now we have the additional modification due to expansion which 
is encoded in both the square root factor as well as the a and /3 coefHcients. There is even a further effect in the 
oscillation frequency due to the renormalization of the frequency represented by We can readily estimate the 

7r|fc| / \^ 2 n> _ZLiM 

values of a and /? for both large and small for large we find Qfe ~ e ( jfy j and Pk ~ ^^p^^ '' 

where the mass index has been suppressed. 



V. THE EQUATION OF STATE FOR THE FLAVOUR VACUUM 

The flavour vacuum as we have seen is orthogonal to the vacuum in terms of mass eigenstates in the thermodynamic 
limit even in the presence of expansion. Arguments based on d-particle scattering of stringy matter give support to 
the flavour vacuum as the correct vacuum. By considering the expectation value of the stress-energy tensor we can 
see to what extent it might play the role of dark energy. This calculation, unlike previous ones, will consider a 
non-flat metric which is of course necessary if there is an expanding universe. Our aim is to investigate the theoretical 
possibility of a dark energy contribution and the issues of D-particle capture and recoil, rather than phenomenological 
relevance. Both early and late time expectation values are calculated. 

The generic stress-energy tensor T^,^ for our scalar theory is formally given by 

Tnu = ]:{(l>,n<i),u + (l),w(l),n)-9nv^= (5.1) 



where L is the lagrangian density, g — det (ff^j/) and = ( the are the generic space-time co-ordinates, 
1-1-1 dimensional in our case). There are of course operator ordering ambiguities which we shall address later. In our 
model 




2 V 9a; 



1 ("^V _ !^^!!l!£M^W2 
2\ dri ) 2 ^ 



(5.3) 



(5.4) 



dr] dx dx drj 



We shall consider expectation values of r^,y for the asymptotic early and late times. For reasons of rigour we will 
work in the discrete momentum basis. In this basis (suppressing in and out indices) let us define 

Oj \r],x) = —a>-'e\^ ^ ' + ay e 37 (5_5j 

(ry,x) ^ a^e^i^^-^'^ + ^Wte-^^-^f") . (5.6) 



and 



In terms of these operators 



-'oo- 2-^00 + 2 "° 2~ °° ' 



-'11-2^00 +2 °° 2~ °° ' ^ 



^i=^EE^}'^?('?'^)^/ (^'^) (5-9) 
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where 



-'oo — 



r(2) 



r(3) 




hf (??, x) 



(5.10) 
(5.11) 
(5.12) 



ana z/^-^-, — 





We will now calculate i„,a,/3(0 T°^* |0) 



/ a,/3,m 



It is straightforward to show that 



t(-) 



and 



,( (^, :r) 4)°"* (ry, :r) |0)„,^,,„ = -%,f5+) + 5, 



(5.13) 



(5.14) 



The expressions for f^-^^ are given in the Appendix . There are infinities in the expectation value of the stress tensor 

which have conventionally been removed by rcnormalization. In flat space-time this has been achieved by a suitable 
normal ordering. Certainly if we consider our model conventionally as a quantum field theory in a time dependent 
metric then we would have to renormalise [10]. This would involve state dependent counter terms which, in a covariant 
procedure, can be tensorially ccmstriictcd from the metric tensor. However out model, as wc have seen, is motivated 
from D-particle capture of stringy matter and so its interpretation cannot be completely that of a conventional field 
theory. We will look at both the low and high momenta behaviour of f^-^^ . The high momenta behaviour will indicate 
the necessity of 'conventional' rcnormalization since the flavour vacuum is indistinguishable from the standard massive 
vacuum the while the low momenta behaviour will give the contribution adapted to the D-particle capture approach. 



The results can be summarised as: 



T, 



(l)otit 



00 



|0)a. 



/3,in 



j j=l 




(5.15) 



(5.16) 



m,a,/3\U -too \^la,l3,in~2L^^ W V 



i I ' 



and 



«,/3(0 Tor|0)„,^,,„ =0. 



(5.17) 



(5.18) 
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Hence 



and 



1 



a,j3,in 2 



^,a,/^(0 



p{l)out 
'00 



|0) 



p(2)out 
'00 



|0) 



+ 



Civ) 



I n T^(3)oMt 
2 in,a,f3 \ " | 00 



|0) 



a,(3,in 



3 i=l 

2 .-2 



J I=i J 

2 



r(l)o"* |n\ 

^00 l"/a,/3,m 



c(/?) 



i=l '^i 



(5.19) 



(5.20) 



The expressions in (5.19) and (5.20) can be normal ordered with respect to a suitable vacuum. Unlike elementary 
field theory there is more than one obvious choice. The equation of state will vary depending on this choice and 
can give physical insight into the fluid like properties of the unusual flavour vacuum in the presence of an expanding 
universe; so a normal ordered expectation can be defined as 

in,c.A^\:nr: io)„,^,,„ = .n,a,/3(o|Tirio)„,^,,„ - {^nrm (5.21) 

where j^*) is a suitable 'vacuum' state. We should remark that various types of 'normal' orderings can be accommodated 
through the choice of |^). The usual one in Minkowski space- time reads: 

: 6« (r;, x) (r;, x) : = 6« (,?, x) (r?, x) + , (5.22) 
: (77, x) d!f (r,, x) : = (r;, x) d!f (ry, x) - 5,,. . (5.23) 

In our D-particle-foam model, which motivates our use of the "flavour" vacuum, the physically correct normal ordering 
is dictated by the underlying microscopic physics of the foamy space-time at Planck scales, and the non-trivial 
interactions with it of some of the flavour momentum modes, specifically the low-energy ones. A cloud of D-particles 
induces a metric which is proportional to (c.f. (2.22), (2.23)). In our picture, the vacuum contribution to the 
expectation value of the stress tensor is due to the D-particle recoil-velocity (statistical) Gaussian fluctuations (2.21), 
i.e. we subtract away the a independent contribution, following our earlier treatment of the induced stochastic light- 
cone fluctuations due to the recoil- velocity fluctuations [6]. This is consistent with the recovery of the Minkowski, 
Lorentz-invariant situation (c.f. (2.22)) in the absence of D-particle recoil-velocity effects, i.e. the vanishing of the 
vacuum expectation value of the matter stress-energy tensor in the limit ct^ ^ 0. 
We should recall here that for a relativistic perfect fluid 



T,j.u = -P9u.v + {p + p) v^Vu 
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where is the 4-vclocity of the observer. From (5. 15), (5. 16), (5.17) and (5.7), (5.8), (5.9) we have the form of the 
expectation values of the various components of the stress-energy tensor, and from those we can conclude about the 
equation of state of our boson gas and check whether p = wp with w = — 1, the expected form for cosmological 
constant. In order to obtain finite expressions for p and p we would need some form of momentum cut-off for 



the low k modes. However, on identifying (in a co-moving cosmological frame) p = 



'0 



TiT-- |0)„,,.,„ and 



'0 



rpout . 

^00 • 



and on adopting the a independent subtraction procedure, we find interestingly (c.f.. 



(5.19) and (5.20)) that 



w 



-1 . 



(5.24) 



It is important to notice that (5.24) has been obtained by a subtraction procedure which was adopted independently 
of the particular choice of the expanding universe conformal factor C(r/). This is important, in that it allows the 
back reaction effects of the flavour vacuum onto the background space time to be incorporated self-consistently in our 
approach. 

The determination of p and p, which are vacuum expectation values, involves a tacit averaging over the time scale 
T in our vacuum. We earlier estimated t to be extremely short and so a time average over it is equivalent to putting 
77 = in the expressions. The contribution to the cosmological constant from this vacuum is given by p. 

It must be noted at this staged that the above calculations have been performed in the regime in which the expansion 
rate of the Universe is slow, and the space-time is almost a Minkowskian one. Hence the cosmological-constant like 
equation of state that characterises our Lorentz-violating vacuum pertains only to late eras of the Universe. In 
general, one might encounter situations with a time dependent equation of state and a time varying dark energy (i.e. 
quintessence). 

We have earlier argued that the capture process is likely to be most efficient for stringy matter with small k. In order 
to get an estimate for the momentum cut-off of the capture process, and thus of the flavoured- vacuum contribution to 
the cosmological constant, we can adopt an argument in [24] by considering the scale inherent in the k dependence 
of the single particle momentum distributions in the flavour vacuum 



n, {k, Tj) = a,/3( al ,^ (??, x) a,,fc {rj, x) |0)„_^ 



for L = a,p. It can be shown that [11] 

na,k {t) = 7ft,- 



sm 



2 • 2 (2) 

■7fe,+ sin ojlJ^kV- 



sin^ 29 



■ sm 



(5.25) 



(5.26) 



and 



n0,k {t) = 7fe,- 

(and we also have the identity 7^ _|_ - 
large k 



sin- C/ 2 • 2 (1) sin 26* . 2 / (i) , (2) \ 
7fe,+ sm LvlJj^T] + sm [u'J,^ + iv'J,^ j r] 



(5.27) 



4). The k dependence is dominated by the behaviour of 7/c,-. Now for 



7k,- 



1 ^(1)2 _ ^(2)2 

2 fc2 



4fc4 



m 



(2)4\ 



and so there is a scale (determined by the ratio of the two terms in the above) 



ko 



1 



(1)2 



7,(2)2^ 



(5.28) 



(5.29) 



which is a plausible cut-off scale in k. We note that fco ~ m -I- | where m^^^ = m + \5m, m^^^ = m 
Hence, there is a fall off for Ua^k [t) with k determined by 7^ _ , which is qualitatively similar to the four-dimensional 
fermion case of [24] and to an earlier conjectured cut-off "^^+™-'^ [18] for rui ^ 1712- 

Prom these considerations, we obtain the flavour-vacuum contributions to the cosmological constant by integrating 
over the extremely small time scales of capture, which according to our previous disc;ussion is qualitatively equivalent to 



setting the time T;-depcndence to 77 = 0. After taking appropriately the contimmm limit, ^ j; Y^j ^ ^ I '^^ 



2"7r lo" with the cutoff ko given by (5.29), we can obtain an estimate for the cosmological-constant contribution 
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due to the flavour- vacuum by considering the case of small expansion and expansion rate. Returning to the expression 
for in,a,f3{0 Tqq^ \^)a f) in (5.19) wc reuormalise the vacuum expectation so as to leave the contribution from the 
string capture of the stringy matter. One part of this renormalisation is the C-independent subtraction which leads 
to 



Stt 
and 



k 



2 ^fen (i)2 



(fff - to • (5.31) 



As noted earlier we already can see the D particle contribution to the flavour vacuum giving w = —1. We further 

reuormalise f^^--* to only include the 7™_ dependent terms since it is only such terms which distinguish it from the 
mass vacuum state. Details are given in the Appendix . To leading order in 5m, on averaging over the uncertainty of 
the capture time 77 (which is effectively instantaneous) and owing to (7.20) and (7.21) 



= 0. (5.32) 

In the next to leading order in 6m 

in,a,0{O ^00'* |0)a,/3,m 

~ iSmf + ^ sin^ e{l-C (,o)) (5.33) 

where r]o is the capture time which in our model can be taken to be rjo = 0. For small a we have C (%) ^ 1 — 45^ 
(c.f. (2.22), (2.23)). Hence, finally making the subtraction of the a independent terms we get 







(Jm)^ 



^oo"*: \0)^.,.,^-'-^smHa'. (5.34) 

This shows the contribution to the dark energy from the flavour-changing D-brane- recoil (statistical) fluctuations. 

The reader should notice that, within our effective field theory framework, the (small) fiuctuation parameter ct^ 
should be considered as phenomenological. In order to determine its real order of magnitude, one needs detailed 
microscopic stringy models of such D-particle foam, along the lines proposed in [4] , in the context of brane- world 
cosmologies (c.f. fig. 1). This issue is still far from being complete, given that realistic three-space-dimensional brane- 
world models of D-particle foam require appropriate compactification, consistent with phenomenologically realistic 
supersymmetry breaking scenarios. We postpone a discussion of such issues for future work. 

A related comment concerns the asymptotic in time nature of the fiavour-oscillation-induced de Sitter vacuum 
energy (5.34). The reader should recall that this result was derived on the assumption that, asymptotically in 
cosmic time, one has a uniform constant density of defects, and the space time approaches the flat, but not standard 
Minkowski (i.e. CT^-dependent), limit (2.22). This will cause eternal acceleration, as the situation is equivalent to that 
of a cosmological constant, and perturbative Scattering matrix and asymptotic states will not be well defined. This 
may have consequences on a quantum field theory of matter defined in such space times, especially from the point of 
view of an i/?-defined CPT operator, a situation which might arise in such a case due to decoherence [25] associated 
with the existence of future horizons (cosmologically induced intrinsic CPT violation [26]). 

On the other hand, in cases where the asymptotic density of defects, crossing our brane world, vanishes, it is 
evident that the flavour- changing induced vacuum energy (5.34) (and also pressure, in view of (5.24) will also vanish 
asymptotically. In such a case, the situation will resemble the one illustrated in the example of the space time (3.4) 
in section III, where the de Sitter phase characterises only a certain era of the Universe, relaxing asymptotically to 
a standard flat Minkowski space time. In the latter case, the perturbative Scattering matrix and asymptotic states 
are then well defined and the associated quantum field theory of matter fields on the brane world is CPT invariant, 
at least from a cosmological view point. 
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VI. CONCLUSIONS AND OUTLOOK 



We have discussed here an approach to the flavour mixing problem and its consequences to the vacuum energy in 
a toy model for bosons, taking into account the effects of the Universe expansion. 

We have attempted to argue that the flavour (Fock space) vacuum introduced by Blasone and Vitiello in order to 
describe canonical quantization of theories with mixing characterises certain string models entailing Lorentz violation, 
in particular D-particle foam models. In such models, the interaction of certain (electrically neutral) string excitations 
with the D-particle defects cause back reaction effects onto the space time, due to the capture of the string state by the 
D-particle and subsequent recoil of the space-time defect. When a statistically significant population of such defects is 
present in our brane world, then their effects are global, and may result in an expanding Robertson- Walker Universe, 
depending on the details of the time profile of the D-particle density. We have provided arguments in favour of such 
a behaviour in the context of specific brane models propagating in a bulk space time, punctured by D-particles. 

During the capture/recoil process, flavour may not be conserved, and it is in this sense that the D-foam background 
provide a way for flavour mixing. The recoil of the defects cause Lorentz violation, as a result of the induced 
momentum-dependent Finslcr-typc [14] metric in the process. Such Lorentz violation may be averaged to zero in 
D-foam models, but leaves traceable effects on quantum fluctuations. 

In such an expanding Universe, with flavour mixing, the vacuum felt by the string states is not the normal vacuum 
but the Fock-space flavoured vacuum introduced by Blasone and Vitiello in theories with mixing. However, in contrast 
to their generic considerations, in our approach the introduction of the flavour vacuum is a result of a specific underlying 
microscopic model of space-time background with defects interacting with some but not all matter excitations. 

The details of the D-particle foam played an important role in determining the characteristic time scales involved 
in the problem (especially the (small) duration of the capture of the string state by the defect), and as a consequence 
the pertinent time intervals over which one has to average in order to obtain the contribution to the vacuum energy. 

We have calculated the effects of the Universe expansion in the flavour- vacuum formalism for a toy model of bosons. 
By taking into account the details of our specific model involving D-particle defects, we have been able of determining 
the physically relevant subtraction scheme in the calculation of the fiavour- vacuum expectation value of the matter 
stress-energy tensor. In this way we have obtained an equation of state of cosmological-constant type for (late) eras 
of the Universe, where the expansion rate is small. This is consistent with the constant contributions to the vacuum 
energy obtained by averaging over the relevant (small) time scales in the problem. Although our considerations are 
toy, as referring to two-dimensional target space time, and hence one cannot make direct contact with dark-energy 
observations, nevertheless, the order of magnitude of this non perturbative contribution (5.34) to the cosmological 
constant can be made phenomenologically consistent with the small value suggested by observations today, provided 
the order of magnitude of the recoil- velocity (statistical) fluctuations cr^, which it depends upon, is sufficiently small. 
Moreover, the contribution vanishes if the mass difference or the mixing angle vanish, as expected in the generic 
approach of Blasone and Vitiello. 

We should extend the above analysis to include fermions, and also we should consider higher-dimensional isotropic 
Universe sitiiations. Moreover, despite the motivations from string theory, our considerations have been restricted 
so far in the context of effective low-energy field theories. It would be interesting to extend these results to pure 
stringy considerations and attempt to estimate microscopically the order of the recoil velocity fluctuations, ct^, in 
realistic models. Such a task is non trivial, as it requires first an understanding of the details of the compactification 
procedure, consistent with supersymmetry breaking in target space. Flux compactification has been briefly used in 
our stringy model in the arguments in favour of the identification of the Liouville mode of our supercritical recoil 
string model with a function of target time. It is our belief that such compactiflcations provide a way of resolving 
the problem of the smallness of the dark energy of the Universe today consistently with phenomenologically realistic 
values of supersymmetry-breaking mass scales. We hope to start tackling such issues in the near future. 

A flnal comment we would like to make concerns the global effect of the induced de Sitter space time on photons. 
According to the analysis in [27], in four-dimensional de Sitter space times, vacuum polarization of charged scalar 
fields, such as charged Higgs in supersymmetric theories, implies mass for the long wavelength photon modes that 
cross the horizon of the (locally) de Sitter space time. Indeed, let k denote the momentum of such modes, 
with c\k\ H initially at some (conformal) time r]o = —H~^, but whose actual (physical) magnitude is diminished 
significantly due to the Universe expansion at a later time rj, after a long period of inflation, \k\/a{r]), with a{r]) 
the scale factor such that a{t) ^ c\k\ ^ H . Then the induced mass m-y of the photon is, to leading order in the 
appropriate expansions [27]: 



where a is the fine structure constant. In our case (5.34), the flavour-vacuum-induced Hubble parameter is extremely 




(6.1) 
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small (at late eras of the Universe), of order \ (Sin)\siir'' OVa^ , and thus only extremely long-wavelength (far infrared) 
photon modes would acquire a (very small) mass via the vacuum-polarization mechanism of [27]. However, in the early 
Universe, our statistical D-particle fluctuations could be much larger, according to the models of D-foam discussed 
above, and thus one does not exclude the possibility of significant photon masses (in supersymmetric theories, where 
charged scalar fields occur naturally) via this mechanism. As discussed in [27], this may produce weak, but of 
cosmological-scale relevance, seed magnetic fields, which could be amplified to produce the micro-Gauss magnetic 
fields observed in galaxies and galactic clusters. It would be interesting to examine such scenarios in the context of 
our flavour-vacuum induced de Sitter space time, as such studies may provide indirect tests of (and impose constraints 
on) space time foam models from novel view points. We hope to tackle such issues in a forthcoming publication. 
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VII. APPENDIX 

In this Appendix we calculate the quantities fi^^, appearing in the appropriately normal-ordered (subtracted) 
expressions (5.30), (5.31) for the boson stress-energy tensor flavour-vacuum expectation value, for slow Universe- 
expansion rate and small momenta k (with respect to the boson mass). 

It is straightforward to demonstrate that 
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We note the general identity 7^ = 4 + 7^ _ and that, for small k, 

Sm 2Am^ m 



7j^ _ ~ ; 

\/m(i)m(2) 4A2 (m(i)m(2))- 



(7.7) 



where A = -'-'j^, -B = ^ and Am^ = m^^-*^ — m^^-*^. In terms of the D-particle model A = 1 — 4(t^ when cr is small. 
The in flavour vacuum differs from the mass eigenstate vacuum through terms proportional to 7^"_ and so the effect 
of the D-particlc capture which we regard as responsible for the mixing in this model is present through such terms 
in f^'^'' and f^"''. Consequently we write for s = 1,2 



and 
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with [1)2 = 1, [2)2 = 2 and [Sjj = 1. The leading contributions to the renormalised : : and : f^^'' : are 
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Now on noting, for real y, that T {iy) T {-iy) = ^ and argT (zy) ^ y\ogy - y - \ + O {j^ as y ^ 00 we can 

show that for small |1 — A| and p that 
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This leads to 
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To lowest order in 5m 



. f(+) . 



COS ^(-^-'-)-(A-i)-^ 

\2p(fc2+m([= + il2)2) 2 



_exp(-lf(Ama-'>J + 53j^))J 



> . 



. . f( + ) . 



(7.18) 



(7.19) 



(7.20) 



. f(-). __..(-). 



and (c.f. (3.34)) 
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which we use in (5.34) to estimate the leading contribution to the dark energy coming from the flavour-changing 
D-brane-reoil fluctuations in our model. 
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